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Effect of in-plane line defects on field-tuned superconductor-insulator 
transition behavior in homogeneous thin film 
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Field-tuned superconductor-insulator transition (FSIT) behavior in 2D isotropic and ho- 
mogeneous thin films is usually accompanied by a nonvanishing critical resistance at low T. 
It is shown that, in a 2D film including line defects paralle to each other but with random 
positions perpendicular to them, the (apparent) critical resistance in low T limit vanishes, 
as in the ID quantum superconducting (SC) transition, under a current parallel to the line 
defects. This ID-like critical resistive behavior is more clearly seen in systems with weaker 
point disorder and may be useful in clarifying whether the true origin of FSIT behavior 
in the parent superconductor is the glass fluctuation or the quantum SC fluctuation. As a 
by-product of the present calculation, it is also pointed out that, in 2D films with line-like 
defects with a long but finite correlation length parallel to the lines, a quantum metallic 
behavior intervening the insulating and SC ones appears in the resistivity curves. 
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1. Introduction 

In homogeneously disordered thin superconducting (SC) films, the resistivity shows an in- 
sulating behavior in higher fields, while it decreases upon cooling in lower fields as a precursor 
of a SC transition. At a critical field separating these two regimes, the resistivity approaches a 
finite value upon cooling. It is well understood that this field-tuned superconductor-insulator 
(FSIT) resistive behavior in two-dimensional (2D) impure superconductors is induced by a 
quantum fluctuation enhanced at lower temperatures. However, two different origins leading 
to the FSIT behavior have been proposed so far. 1-4 One is the quantum critical fluctuation 
accompanying the continuous melting transition of the 2D vortex-glass (VG) state at T = 0. 
The other is a consequence of the thermal to quantum crossover behavior in the conventional 
Aslamasov-Larkin fluctuation conductivity cr s jj (j = x, y): The thermal fluctuation enhances 
cr s jj upon cooling, while the quantum fluctuation is insulating 3 and reduces <r s jj upon cool- 
ing. Since, as the magnetic field is lowered, the fluctuation is weakened, and the quantum 
fluctuation is changed into the thermal one, just the SC fluctuation itself may induce the 
FSIT behavior. However, this picture, first argued in Ref.3, was found not to be justified once 
the vortex pinning effect is taken into account at the microscopic level for the homogeneously 
disordered SC films. 5 Nevertheless, it has been argued recently that the FSIT behavior in 
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Fig. 1. (a) Sketch of a SC thin film under a magnetic field parallel to z and with line defects (dashed 
curves) parallel to the y axis, (b) The corresponding figure including line- like defects with a finite 
correlation length along the y axis in place of the line defects in (a) . 



NdCuO is a consequence of the conventional fluctuation conductivities including cr s jj. Un- 
der such situation, it will be valuable to propose an experiment for judging which of the two 
origins of an FSIT behavior plays a dominant role. 

In this paper, effects of line defects on resistivity curves in a homogeneously disordered 
thin SC film are considered at relatively low temperatures. The film sample is assumed to 
be isotropic prior to introducing the line defects running along with one another, and their 
positions perpendicular to their extended direction are not periodic but random. As sketched 
in Fig. 1(a), such line defects may be artificially introduced in real amorphous thin films. It 
is found by examining the SC part of conductivity that the line defects make the scaling 
behavior of resistivity anisotropic : The line defects merely affect the magnitude of the SC 
part of conductivity when the current is perpendicular to the line defects, while they make the 
nonvanishing conductance 1 ' 2 at the (apparent) critical field H c divergent for a current paral- 
lel to the line defects. We find through numerical analysis of resistivity that this anisotropic 
behavior appears more remarkably in the case with a smaller R n (i.e., weaker point disorder), 
where the flat resistivity curve at a critical field for the perpendicular current is vague re- 
flecting that the FSIT behavior originates largely from the quantum SC fluctuation, while the 
vanishing critical resistance for the parallel current is well defined. Thus, the present result 
is relevant to understanding the origin of the FSIT behavior in real systems correctly. As a 
byproduct, we also examine the case such as Fig. 1(b) including line-like defects with a finite 
correlation length and find that, reflecting a ID to 2D crossover upon cooling, the resistivity 
curves in currents parallel to the line-like defects, close to the critical field, decrease rapidly 
at intermediate temperatures, while they become flat at the lowest temperatures even in sys- 
tems with weak quantum fluctuation. Implication of such a nonmonotonic behavior will also 
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be discussed in relation to real observations. 6-9 

In sec. II, analytical expressions of conductivities are derived on the basis of the Ginzburg- 
Landau (GL) model. In sec. Ill, numerical results of the conductivities in sec. II are shown, and 
their contents are discussed. Consequences of the present results are discussed in sec. IV. 

2. Model and analytical calculations 

To explain the analysis used in this paper, let us first start with the 2D Ginzburg-Landau 
action 5ql to be derived microscopically 

f/3 



Sgl = sj d 2 r|/3^C(r)7(Q 2 )kl^(r) + ^ dr 
+ ^(r, r)| 4 + u(r)p(r, r) + /(r)p„(r, r) 



V>*(r, r) / u(Q 2 )V'(r, r) 



(1) 



where p(r, r) = |V(r, r)| 2 , p v (r, r) = [Vx J(r, r)] 2 , J = fg(^*Q^+c.c.), Q = -iV + 27rA/</» 
is the gauge-invariant gradient, and random potential terms leading to pinning effect of vortices 
were assumed to be present. As shown elsewhere, in the absence of line defects lying in the 
plane, the corresponding replicated action rather than eq.(l) itself is derived microscopically. 
Here, 4>{r, r) = Y^, w V^( r ) e~ 1WT , b > 0, s is the film thickness, r is an imaginary time, and 
u) is a Matsubara frequency. Further, when a small electron-hole asymmetry is neglected, 7 
is real and positive so that the dynamics of the pair-field is purely dissipative. The random 
potentials will be assumed to satisfy the following Gaussian ensembles 

u(r)u(r') = w(r - r') = A p 5(r — r) + ^-5(x — x'), 

w(k) = [ d 2 rw{r)e-^ r = A p +'^6(k y ), 
J ?o 

WW) = w^r-r'), 

«;*(k) = / d 2 r W$ (r)e- ik ' r = + —J^8(ky), (2) 
J ?o 

where £0 is the coherence length defined in the ordinary GL region (i.e., in low fields and near 
T c ). Hereafter, the length scale and the pair-field will be treated in dimensionless forms under 
the scale transformations r/r# — ► r, and ip (/3s) 1//2 — > ip, where rn = ^/ (J)o/(2ttH). 

After performing the random average, we encounter a replicated action S n corresponding 
to eq.(l), where 

S n = S% + S;, (3) 



and the free energy is given by 



and 



7 n — 1 

T= lim -, (4) 

n 



Tr^e- 3 ". (5) 
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The first term of eq.(3) is the contribution independent of the pinning effect and, using the 
Landau expansion of the pair-field 

^ Q) (r) = E^WK.pM, (6) 

n,p 

becomes 



So = EE (^o+Toi^DEi^^^i' + ^iHr'Ei^^^) 

a u> L p p 



(7) 



where a is a replica index, 

4 Q) (k, c) = E exp Apfc, - ^) (4°V, ^'))Vo a W " + a/) (8) 

(p± = p± fej,/2), and 

g 1 H = (| 9 9 1 (p,u;)| 2 > = ( /Ul +7 1 |u;|)- 1 . (9) 

Here, u„ jP (r) is the eigenfunction in the n-th LL, and only the lowest (n = 0) and the next 
lowest (n = 1) LLs were kept. Further, N v is the number of field-induced vortices, and the 
dependences of 7„ and fj, n on the LL index n stem from the Q-dependence of the corresponding 
coefficients in eq.(l). 

Further, note that, in eq.(9), [i\ is the bare mass of the n = 1 LL fluctuation. That is, we 
have included the n = 1 LL fluctuation only at the Gaussian level. Since this fluctuation is 
heavy, and fj,i is of order unity in high fields and low temperatures relevant to the quantum 
fluctuation phenomena, the renormalization of n = 1 LL fluctuation due to interactions among 
fluctuations is negligible. For this reason, the Gaussian approximation for the n = 1 LL 
modes is adequate. On the other hand, the lowest LL fluctuation should be treated in a 
fully renormalized form, because the events in equilibrium, such as the vortex lattice melting 
in clean limit and the glass transitions, occur in the lowest LL in the present high field 
approximation. For the present purpose of examining resistive behaviors in the vortex liquid 
region, the one-loop approximation used elsewhere 5 for the lowest LL fluctuation propagator 
Gq(uj) = (\<fo{p, ^)\ 2 ) w iH be sufficient. Then, ^o(^) is given by 

SbM = <|rf*W)| 2 >, 



( 70M+ , 0+so -IM( AP+A ^))-. (10) 
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where h = Co/ r H = ^nH^/cfto. The pinning-induced contribution S™ to the action (3) is given 
by 



" -^EEk k )^ (Q) ( k ' °))*/5 (/3) ( k ' 0)+^(k)(pH(k, 0))*pf(k, 0) 



a,/3 k 



(12) 



where 



p(k, 0) = y*d 2 rp(r, 0), 



— ikr 



<Po(p+, wM>(P-> w) - i^=¥>I(p+, w)<p (p~, w) 



(13) 



p„(k,0) = /iE^kx) 172 ^ 



fc-. 



(14) 



— dz iky j cincl — /c^ "I - . 



Therefore, 5™ takes the form 



A„(k)p<°»(k, 0)p<»(-k, 0) 



where 



o, /3 k 

+ ^AxCkJ^Oc, 0))*pf >(k, 0) + ^A x (k)pf )(-k, 0)(pf } (-k, 0))* 
+ ^(kJpfVk, 0)pf (k, 0) + f MkXpf )(k, 0))*(pf)(-k, 0))* 
+ i^A 01 (k)p( Q) (-k, 0)pf } (k, 0) - i^A 01 (k)(pf } (k, 0))*pf (k, 0) 

- i^A 01 (k)pf } (-k, 0)pf (k, 0) +i A_A 01 (k)4 Q) (-k, 0)(pf } (-k, 0))* 



Ao(k) = expf-M 

Ai(k) = expf—) (A p + /i 2 A<*>(fc 2 -2) 2 ) 



(A p + tf A (* } fc 4 ) + ^(A, + h^t 4 )^) 



(15) 

(16) 
(17) 



+ ^(A, + ^ 2 A^(fc 2 -2) 2 )^) 



(18) 



5/16 



J. Phys. Soc. Jpn. 



Full Paper 




Aoi(k) 



exp 



2 



(A p + h 2 A^k 2 (k 2 -2)) 
(Aj + h 2 Af ] k 2 (k 2 -2))6(ky 



(19) 



In the absence of line defects, i.e., when A/ = A^ = 0, the above expressions can be 
regarded as being derived microscopically. Describing line defects introduced artifitially from 
the microscopic standpoint is not easy, and they were incorporated in a phenomenological but 
conventional manner. 
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We examine hereafter the superconducting part cr s ,jj of conductivity <Tjj (j = x, y), ac- 
cording to the Kubo formula, in the form of a current-current correlation function by invoking 
the high field approximation where the expression of the supercurrent density is linear both 
in the n = and 1 LL modes of the pair-field. In the framework consistent with the one-loop 

approximation for the lowest LL fluctuation, a s ,xx takes the form 

2 



P,P',U! 



(20) 



(flip, u + ^)<Pq(p', v + ^■)){<fi(p', u)<Pq(p, w)> 

while the corresponding expression of <J s ,yy is given by replacing the minus sign prior to 
the second term of eq.(20) by a plus sign. Here, Rq = 2irh/Ae 2 is the quantum resistance. 
Note that, because of line defects introduced along the y-axis, the second term of eq.(20) is 
nonvanishing. 

In the presence of pinning disorder, a s> ^ can be seen as the sum of the so-called vortex 
flow term <7fl )W , independent of the vortex-glass fluctuation, and the divergent term cr vg)Wi 
reflecting a growth of the vortex-glass fluctuation. A method of formulating o" vgj w was devel- 
oped in Ref.10 in the case with only point defects and with no quantum SC fluctuation. As far 
as the lowest order contribution in is concerned, where £ vg is the glass correlation length 
normalized by rjj , ^vgjj obtained in the lowest order in the pinning strengths consist of the 
sixteen diagrams in Fig. 2. The second term of eq.(20) is given by the sum of diagrams 1)' to 
8)'. In the expressions of these diagrams, the vortex-glass correlation function G vg (k ; u, lo + SI) 
inducing the £ V g-dependences is included. In the present case with line defects perpendicular 
to the field, G vg has been studied elsewhere 10 and, in the so-called ladder approximation, 
becomes 

(So(0)) 2 



£? vg (k; uj,(jj + CI) 



r 9 + E^y^k 2 + r 9 So&W (M + \uj + SI 

t-2 m ln\\2 



where c y ~ 1/2, 



l + ^ g E c^ + 4 7o a (0)|a;|+47oao(0)|a; + fi| 

H=x,y 



c x ^- Ap , (22) 

27T 
h 



2 A p + A, 




and 



As shown in Ref.10, this glass fluctuation in 3D systems is 2D-like if the point disorder is 
absent. Reflecting this feature, the coefficient c x vanishes in A p , A p ^ — > limit. 
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Diagrams contributing to the r.h.s. of eq.(20) are described in Fig. 2, where each open 
rectangle corresponds to G vg . Eight pairs of diagrams shown in Fig. 2 take the form 



l) + 2) 
l)' + 2)' 

3) +4) 
3)' + 4)' 
3)' + 4)' 

5) + 6) 
5)' + 6)' 

7) + 8) 
7)' + 8)' 



k + k'_ 



, 2 



V(k)V$(k') 







k Jv 2 

k+k'_ 

k J\u 



k+k>+ V(k)V*(k') 



d\n 

d 



F(k)y$(k') 
y(k)T4(k') 
F(k)y$(k') 
y(k)T4(k') 
y(k)T4(k') 
F(k)T4(k') 



d\n 

d 



k + k' + _ 

k + k' + 

k+k'_ 
'kJk' 2 

k + k' + 

■ ./k' 2 

k+k'_ 
k Jin' 2 



9|n 



9|n 
5 



9|n 



9|n 



9|n 
o 



k Jk 



, 2 



d\n 

d 



a\n 



/T 1 ]T(g (u,)) 2 g vg (k - k'; w , w + n), 

UJ 

/T 1 o M0o(w + ^)a v g(k - k'; u, uj + fi), 

UJ 

/T 1 ^(C?o(^ + ft)) 2 a vg (k - k';w,w + 0), 

UJ 

/T 1 g (^)^o(^ + ft)£v g (k - k'; w, w + ft). 

UJ 

I 3 ' 1 S £oM£o(^ + ^)a vg (k - k'; w, w + fi). 

UJ 

/r 1 Yj ^)Go(oo + n)g vg (k - k' ; uj, u + n), 

UJ 

/T 1 Y,(So(^)) 2 G vg (k - k'; w , w + 17), 

UJ 

I 3 ' 1 S &H0o(w + ^)£?vg(k - k'; w, w + O), 
/T 1 ^(C?o(^ + ^)) 2 a v g(k - k';w,w + 0), 



where we have used the fact that the wave number k — k' carried by the glass fluctuation is 
small so that 



and 



Ai(k)A (k') - Aoi(k)A i(k') ~ 4h 2 V(k)V*(k'), 



V*(k) = 



exp 



A; 2 \ 2tt 



v~h Al exp 



S(ky). 



(24) 



(25) 



(26) 



Further, the frequency dependence of the high energy fluctuation in the next lowest LL was 
neglected so that Q\ = (|v?i| 2 ) is the constant /i 1 " 1 Cbecause we are interested in the time scales 
of the glass fluctuation with the lowest energy in the present situation. 

Below, we focus on the real part of the conductivities. Then, Rq a vgtXX is given by 

!) + ••• + 8)- (!)' + ••• +8)' ) 



jf ^ k y k' y V(k)V$ (k') (-^| ) P' 1 + + ^))^vg(k - k'; u, uj + Q) 

/r 1 Y(G (uj)+go(Lu + n)) 2 I g vg (k';u,uj + n), (27) 



h 2 J^A p A^k 2 e 



8/16 



J. Phys. Soc. Jpn. Full Paper 

while the corresponding expression of Rq cr vg: yy is 
l) + --- + 8) + ( 1/ + --- + 8)') 

= h 2 JJ^ k x k' x V (k)v*Qi')(--^Jp- 1 52(g (u) + g (u; + n)) 2 g vg (k - + n) 

~ /i 2 A p Af e" fc2 (-^| ) r 1 J>oM + ^ + ^)) 2 £ ^g(k'; ^ + 0) 
x / e -^ 2 / 2 g vg (k / ;L I j,L I j + ^) 

+ hj^ A,Ap ) fc^- fc ^-^/3- 1 ^(g H + g (a; + n)) 2 ^ / av g (4;^^ + ^)- (28) 

3. Numerical analysis of resistivity curves 

To discuss the resistive behaviors in the quantum regime in details, we examine the re- 
sulting resistivity curves numerically. Before performing this, the expressions of conductivities 
need to be arranged in a useful form. 

First, strictly speaking, the pinning strengths such as A p and A p ^ carry the k or k' 
dependences. However, they are not obtained in a closed form. In particular, it is not easy 
to know the momentum dependence of A; and A[*^. For these reasons, we write, e.g., the 
coefficient J k h 2 A p A p ^k 2 exp(-fc 2 ) in eq.(27) as c p (A p ) 2 /(r^s 2 ). Note that A p /(^s) is a 
dimensionless quantity. Since we expect the magnitude of A p ^ to be much smaller than that 
of A p , the coefficient c p is assumed to be smaller than unity. NextCarranging the frequecy 
summation using the replacement Yl a <ui<b y ^( aJ ) — ¥ S^> a ~~ Ylu>b y ^( u; ) an d performing 
the k'-integral over the glass correlation function, the expression of cr vgi xx becomes 



£ ^ [ g0H(60( ; ) + g ° (0)) ln(l + *2 



l + k 2 c (c x c y ) l / 2 e v g + 2^| 70 g (0)^ 



(29) 



l + 2M 7 o<?o(0)£4 g 

where k c is a cut-off of order unity of the |k|-integral. Further, the prefactors 
J kx h 3 / 2 (A p A^ ] + AiAf ) )fc 2 exp(-fc 2 ) and j kx hA^ k 2 exp(-k 2 ) will also be replaced 
with 2h~ 1 / 2 A p AiC\ p /(rf]S 2 ) and h~ 1 (Ai) 2 c\/(rj I s 2 ), respectively, by introducing unknown 
constants q p and q measuring nonlocalities of the pinning potentials due to line defects. 
Then, the conductivity cr vgiTO takes the form 

sR Q a vg , yy = sR Q * vg , xx (l + 2^h-^ +^ g ^^(l + ^/f ) V V 1 (a (0)) 2 
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Fig. 3. Examples of resistivity curves computed in terms of eqs.(29) and (30), R xx ((a) and (c)) and 
R yy ((b) and (d)). The curves in (a) and (b) are those when R n — 0.39-Rq, and the values of the 
reduced field h are h = 0.69 (bottom), 0.74, 0.77, 0.79, 0.81, 0.84, 0.88 (top), while the curves in 
(c) and (d) are those when R n = 0.52i? Q , and the h values 0.60 (bottom), 0.65, 0.68, 0.70, 0.72, 
0.75, 0.80 (top) were used. Each curve at H = H c is indicated by a solid curve in each figure. 
Other parameters used here are A//(£qs) = 3.0 x 10~ 3 , c p = 0.2, q p = 0.2 and c\ = 5.0. 



Y — 



g M(goM + go(o)) _ (goH)2(1 + 2 , a; | 7ogo(0)4) -i/2' 



(30) 



2 

We guess that, compared with the point defects, real line defects will have a longer correlation 
range in the x direction, because they, if existing at low densities, have a strong impact. 
Hence, we expect c\ to be somewhat larger compared with c p . When focusing on metallic 
superconducting films to be described in dirty limit, A p is identified with b p defined in Ref.2, 
while the dimensionless strength Aj/(£qs) of the artificial line defects is unknown and will be 
assumed to be a constant independent of the magnetic field. 

Besides these expressions of the glass fluctuation terms of conductivities, we have included 
the contributions insensitive to the defects to the total conductivities, the quasiparticle term 
a n and the ordinary SC fluctuation term or both of which were defined in Ref.12. In eq.(10) 
giving the SC fluctuation propagator G(u>) selfconsistently, a fluctuation renormalization of the 
bare pinning strengths A w (w = p, /) will be incorporated by replacing them by A w /(1 + v p ), 
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where 

^ = |^£(So(M)) 2 - (31) 
Consistently with this treatment, the glass correlation length is give by 

^^-^h( A ' + A '^) + ^ rm " (32) 

where the second (~ T 1 / 2 ) term, which should appear in the vicinity of the quantum critical 
point of the glass transition, was incorporated. 2 

Using the above expressions of each term of the conductivities, together with the expres- 
sions of /i(0), b, A p /(r 2 js), and 70 in dirty limit, 12 we have numerically examined the resulting 
resistance v.s. temperature curves in order to clarify the significant roles of line defects. Below, 
R n denotes the resistance in the normal state of a superconducting film with point defects 
but no line defects. We choose the parameter values for the dimensionless strength of line 
defects A//(£ 2 s) = 1.0 x 10~ 3 , 3.0 x 10~ 3 , or 4.5 x 10~ 3 , the momentum cut-off k c = 1.0, 
R u /Rq = 0.52, 0.47, or 0.39, c p = 0.2, q p = 0.2, and q = 5.0. For the reason mentioned 
earlier, a much smaller value than q is assumed for c p . 

Resistivity curves following from eqs.(29) and (30) are shown in Fig. 3. As expected, the 
conductivity a xx for a current perpendicular to the line defects includes no explicit effects of 
line defects, as expected from the mean field picture that the vortex flow motion in this case is 
unaffected by the presence of line defects. Actually, the resistivity curves in a narrow window 
of the magnetic field around an estimated critical field H = H c show a flat behavior insensitive 
to the temperature at low enough temperatures, although, as already seen elsewhere, 12 the 
flat behavior near H c is not visible in the case with low R n where the quantum fluctuation is 
weaker. We note that, for R n = 0.52R Q , H c = O.7O0o/(2vr£ 2 ), while H c = 0.79(/> /(2< 2 ) for 
Rn = 0.39Rq. On the other hand, the response a yy to a current parallel to the line defects has 
a term proportional to the square of the strength of line defects, and this term proportional to 
q is dominant near the critical point at T = where £ vg diverges. In fact, in situations where 
the glass fluctuation plays a significant role, the glass fluctuation terms oVgjj of conductivities 
show the following behaviors 

R Q sa^ yy ^^ g Y{T^), (33) 

where z vg is the dynamical critical exponent of the glass fluctuation and, in the present mean 
field (i.e., Gaussian) approximation, is four. In particular, just on the critical field defined at 
low enough temperatures, i.e., H = H c , we have the relation £ vg ~ T~ l l Zv % (see also eq.(32)) 
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when T — ► 0. Then, the critical value of cr vg jj becomes 



Rq & Vg, XX 



Rq OVg, yy 



(H = H c ) ~ (nonuniversal) constant 



(34) 



(35) 



This divergent a V g !yy (H = H c ), indicative of a vanishing critical resistance, corresponds 
to that 13 expected in the virtual ID superconductor-insulator quantum transition where the 
critical conductance is infinite. As seen in the figures, this ID behavior in a yy is more visible 
for a high density of line defects and/or for a smaller R n /RQ value (i.e., weaker quantum 
fluctuation). As the figures suggest, the critical field seems to be well estimated from an 
onset field of the saturated behavior of R XX (T) curves, while in the case with low R n /RQ 
the critical field is not well estimated from the corresponding R XX (T) curves in which the 
saturated (flat) behavior is vague and not clearly seen. That is, when the flat behavior in a 
vanishingly narrow field range is clearly visible in R XX (T) curves, such a flat behavior should 
be a direct consequence of the quantum vortex-glass fluctuation, while the position of critical 
field in systems with low R n /RQ is more clearly seen in R yy data by artificially introducing 
line defects in the y-direction. 

4. Discussion 

As an application of the present result, it is interesting to imagine the situation illustrated 
in Fig. 1(b) where, instead of the line defects, line like defects with a long but finite correlation 
length l y of its orientation parallel to the y-axis are present in a SC film. As in the case with 
line defects parallel to the y-axis, there are no essential changes in a xx , and effects of the finite 
l y appear mostly in cr vg , yy . This finite l y is naturally incorporated by multiplying the last term 
ofiw(r — r') (w<s>(r — r')), proportional to A; (A^), by exp( — (y — y') 2 /{2 ly) ). Then, as far as 
A;y / 27r//i <C Ap, effects of the finite l y on Qq{oj) and £ vg are unimportant, and the essential 
change can be incorporated by replacing the last line of eq.(30) with 



where l y = l y /rn- 

From this expression, the following qualitative behaviors are expected on R yy . At higher 
temperatures where £ vg <C l y , the ID behavior, eq.(33) or (35), of cr yy {H ~ H c ) is valid, and 
the R yy (T) curve decreases upon cooling as in Fig.4(b) in fields around and below H c . In 
contrast, at lower temperatures where £ vg 3> l y close to or below H c , the growth of cr vgjyy due 
to the line-like defects upon cooling saturates in some field range around H c , reflecting the 
ordinary 2D behavior. In fact, as the R yy data in Fig.4(c) following from eq.(30) with eq.(36) 
show, such resistivity curves around H = H c ~ 0.74(/>o/(27r£g) that are flat at the lowest 



x 




d r goH(goH+g (o)) 

d\(j\[ 2(l + g 2 ) 1 /2 



(SoM) 2 



(36) 



(l + g 2 + 2M 7o So(0)e)i/2j 
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0.10 0.20 0.30 0.40 



Fig. 4. Computed (a) R xx curves, (b) R yy (T) curves in the case with persistent line defects parallel 
to y, and (c) R yy (T) ones in the case with line-like defects with l y = 0.75 (i.e., in the situation 
sketched in Fig.l (b)) at the helds h = 0.73, 0.733, 0.735, 0.737, 0.739, 0.74, 0.741, 0.742, and 
0.744. The parameter values used here are R n = 0A7Rq, A ; /(^s) = 3.0 x 10~ 3 , c p = 0.2, q p = 0.2 
and c\ = 5.0. 

temperatures but decrease upon cooling at higher temperatures become more remarkable 
when, as in Fig. 1(b), line- like defects with a finite l y are included. This feature arising by 
assuming the presence of highly anisotropic defects may be relevant to the data 6,7 suggesting 
the presence of an intermediate "bose-metal" phase. 

Further, when the strength of line-like defects is stronger, reentrant resistivity curves, 
showing a decrease in intermediate temperatures but an upturn in lower temperatures, are 
more remarkable, just as in the curve at h > 0.76 in Fig. 5(c), because stronger line defects 
may induce a sharp decrease of resistivity even in H > H c (i.e., h > 0.746). Such an observed 
reentrant resistivity curve, seen in earlier data of nominally amorphous materials, 8 is often 
regarded as implying that the sample is a granular film. 14 Although, roughly speaking, the 
present picture on the reentrant behavior is similar to assuming a granular structure in the 
film sample in the sense that a mesoscopic structure is invoked, we believe that the present 
picture is realistic in the sense that a local anisotropy of inhomogenuity is incorporated. 

In addition, we note that, in the case corresponding to Fig. 5(c), the critical resistance 
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Fig. 5. Computed resistivity curves which are the same as those in Fig. 4 except the use of A;/(£gs) = 
4.5 x 1(T 3 . h = 0.737, 0.741, 0.744, 0.746, 0.748, 0.749, 0.75, 0.751, 0.753, 0.755 and 0.758 respec- 
tively. 




Fig. 6. Sketch of a SC film including line-like defects with randomness in their directions. 



value, which is ~ 0.5R n in Fig. 5(c), is highly suppressed by the presence of l y of order unity. 
As already pointed out in Ref.12, the apparent value of the critical sheet resistance is close to 
Rn in most cases. However, some data show a significantly low critical resistance especially 
in systems with a larger R n /RQ. 9 The present result will give one of possible explanations of 
such a significantly low critical resistance estimated experimentally. 

Although an appearance of the "bose-metal" behavior or a reentrant one has been demon- 



14/16 



J. Phys. Soc. Jpn. Full Paper 

strated above for the case with line-like defects directed along a unique direction (see Fig. 1(b)), 
similar behaviors should also be seen in the case, as in Fig. 6, where the correlated direction 
of line-like defects is random, because the regions with line-like defects parallel to the applied 
current have the largest weight on the total conductivity and thus, play dominant roles in the 
in-plane resistivity curves for any current direction. Therefore, if the film sample is globally 
isotropic but has such a highly anisotropic inhomogeneous structure over a finite length scale 
longer than the averaged vortex spacing, we argue that unexpected resistivity curves showing 
the "bose-metal" behavior or a reentrant one may appear, reflecting the ID to 2D dimensional 
crossover induced by the line-like defects correlated over such a finite length scale. 

The computation in this work has been done using the facilities of the Supercomputer 
Center, Institute for Solid State Physics, University of Tokyo. 
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